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Abstract. First-principles, density-functional-based electronic structure calculations 
are employed to study the changes in the electronic properties of ZnCyNis 
and MgC y Ni 3 using the Korringa-Kohn-Rostoker coherent-potential approximation 
method in the atomic sphere approximation (KKR-ASA CPA). As a function of 
decreasing C at%, we find a steady decrease in the lattice constant and bulk modulus in 
either alloys. However, the pressure derivative of the bulk modulus displays an opposite 
trend. Following the Debye model, which relates the pressure derivative of the bulk 
modulus with the average phonon frequency of the crystal, it can thus be argued that 
ZnCNi 3 and its disordered alloys posses a different phonon spectra in comparison to 
its MgCNiz counterparts. This is further justified by the marked similarity we find in 
the electronic structure properties such as the variation in the density of states and the 
Hopfield parameters calculated for these alloys. The effects on the equation of state 
parameters and the density of states at the Fermi energy, for partial replacement of 
Mg by Zn are also discussed. 
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1. Introduction 

In spite of being iso-structural and iso-valent to the cubic perovskite 8K superconductor 
MgCNi 3 [T], ZnCNi 3 remains in the normal metal state down to 2K [2]. The specific 
heat measurements indicate that the absence of superconductivity in ZnCNi 3 may 
be due to a substantial decrease in the density of states at the Fermi energy N(Ep) 
resulting from its relatively low unit cell volume in comparison with MgCNi 3 [2j. 
However, electronic structure calculations show that the decrease in N(Ep) is not sizable 
enough to make ZnCNi 3 non-superconducting [H]. For both MgCNi 3 |U [3 El Ej and 
ZnCNi 3 [H] the density of states spectra display similar characteristics, particularly in 
the distribution of electronic states near the Fermi energy Ep. The electronic states at 
E F are dominated by Ni 3d states with a little admixture of C 2p states. There exists a 
strong van Hove singularity-like feature just below Ep, which is primarily derived from 
the Ni 3d bands. 

To account for the lack of superconductivity in ZnCNi 3 , the density- functional 
based calculations emphasize that the material subjected to the specific heat 
measurements may be non-stoichiometric in the C sub-lattice PJ. This would then make 
it similar to the a— phase of MgCNi 3 , which has a low unit cell volume and remains 
non- superconducting [E]. It has been shown earlier that exact C content in MgC y Ni 3 
depends on the nature of synthesis and other experimental conditions [Tl IH1 IH1 ITTH fTT| IT^ . 
According to Johannes and Pickett p], the arguments that favor non-stoichiometry are 
the following: (i) Total energy minimization en-route to equilibrium lattice constant 
within the local-density approximation (LDA) finds an overestimated value for ZnCNi 3 
in comparison with the experimental values. In general, overestimation is not so common 
in LDA. Meanwhile, when one uses similar technique for MgCNi 3 , the calculations find 
a slightly underestimated value which is consistent within the limitations of the density- 
functional theory lUHHUHj- (ii) The authors also find N(Ep) in MgCNi 3 estimated 
as 13.6 states/Ry atom, while for ZnCNi 3 , under similar approximations, it was found 
to be 11.01 states/Ry atom. Note that it has been shown both experimentally as 
well as from first-principles calculations that a decrease in the lattice constant or a 
decrease in the C occupancy would lead to a decrease in N(E F ) (THj. (iii) A decrease 
in the unit cell dimensions can induce phonon hardening. This is well supported by 
the experiments which find the Debye temperature approximately 1.6 times higher for 
ZnCNi 3 in comparison to MgCNi 3 [2\. 

Earlier synthesis of ZnC y Ni 3 [13 [THJ [TJj finds the lattice constant to be 6.899 
a.u., for which the occupancy in the C sub-lattice was just 70%. The authors have 
employed similar preparation technique for MgCNi 3 [T^] and have found that the 
C occupancy ranges between 0.5-1.25, which is consistent with the recent reports 
PJ El EH EH CH 112 HE|- Lattice constant for ZnCNi 3 , as high as 7.126 a.u. has also 
been reported elsewhere [HI EH], which then becomes consistent with the recent total 
energy minimized value using density-functional based methods. Hence, it seems that 
ZnCNi 3 which was subjected to specific heat experiments [2 J may indeed suffer from 
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non-stoichiometry. 

To understand and compare the effects of C stoichiometry on the structural and 
electronic properties of MgC y Ni 3 and ZnC y Ni 3 , we carry out a detail study using 
the Korringa-Kohn-Rostoker (KKR) Green's function method (2TJ [22] formulated in 
the atomic sphere approximation (ASA) [23J. For disorder, we employ the coherent- 
potential approximation (CPA) [2lj. Characterization of MgC y Ni 3 and ZnC y Ni 3 with 
0.85 < y < 1.00 mainly involves the changes in the equation of state parameters viz., the 
equilibrium lattice constant, bulk modulus and its pressure derivative. The electronic 
structure is studied with the help of total and sub-lattice resolved density of states. 
The propensity of magnetism in these materials is studied with the help of fixed-spin 
moment method [2H] in conjunction with the Landau theory of phase transition [26j. 
The Hopfield parameter g which generally maps the local "chemical" property of an 
atom in a crystal is also calculated as suggested by Skriver and Mertig and its 
variation as a function of lattice constant has also been studied. In general, we find 
that both MgCNi 3 and ZnCNi 3 display very similar electronic structure. Evidences 
point that the non-superconducting nature of ZnCNi 3 may be related to the crystal 
structure characteristics, namely phonon spectra. 

2. Computational details 

The ground state properties of MgC y Ni 3 and ZnC y Ni 3 are calculated using the KKR- 
ASA-CPA method of alloy theory. For improving alloy energetics, the ASA is corrected 
by the use of both the muffin-tin correction for the Madelung energy [2B] and the 
multi-pole moment correction to the Madelung potential and energy [221 f3Ej]. These 
corrections have brought significant improvement in the accuracy of the total energy by 
taking into account the non-spherical part of polarization effects [SI]. The partial waves 
in the KKR-ASA calculations are expanded up to l max = 3 inside atomic spheres, 
although the multi-pole moments of the electron density have been determined up 
to l^ ax = 6, which is used for the multi-pole moment correction to the Madelung 
energy. In general, the exchange-correlation effects are taken into consideration via 
the local-density approximation with Perdew and Wang parametrization [22], although 
a comparison in the equation of state parameters has been made in this work with the 
generalized gradient approximation (GGA) [33j. The core states have been recalculated 
after each iteration. The calculations are partially scalar-relativistic in the sense that 
although the wave functions are non-relativistic, first order perturbation corrections to 
the energy eigenvalues due to the Darwin and the mass- velocity terms are included. The 
atomic sphere radii of Mg (Zn), C and Ni were kept as 1.404, 0.747, and 0.849 of the 
Wigner- Seitz radius, respectively. The vacancies in the C sub-lattice are modeled with 
the help of empty spheres, and their radius is kept same as that of C itself. The overlap 
volume resulting from the blow up of the atomic spheres was less than 15%, which is 
legitimate within the accuracy of the approximation [34J. 

The electron-phonon coupling parameter A can be expressed as g/M (a> 2 ), where 
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rj is the Hopfield parameter, expressed as the product of N(E F ) and the mean square 
electron-ion matrix element (J 2 ), with M and (lo 2 ) being the ionic mass and average 
phonon frequency [35J. However, one may note that the above decomposition of the 
problem into electronic and phonon contributions is only approximate since in principle 
(co 2 ) is also determined by the electronic states. It follows that the Hopfield parameter is 
the most simple basic quantity which one may obtain from first-principles as suggested 
by Gaspari and Gyorffy The latter assumes a rigid muffin-tin approximation 

(RMTA) in which the potential enclosed by a sphere rigidly moves with the ion and the 
change in the crystal potential, caused by the displacement, is given by the potential 
gradient. Within the RMTA the spherically averaged part of the Hopfield parameter 
may be calculated as, 

V ° = N(E7) ? (/ + l)Ml ^ (2/ + l)(2/ + 3) (1) 

where N(E F ) is the total density of state per spin at the Fermi energy, and N t the I th 
partial density of state calculated at the Fermi energy E F , on the site considered. The 
term Mij + \ is the electron-phonon matrix element given as |27j . 

M M+1 = / r'Rt—R^dr (2) 
Jo dr 

which are obtained from the gradient of the potential and the radial solutions R[ and 
Ri+i of the Schrodinger equation evaluated at E F . The special form of the EqsH]and[2] 
stems from the ASA in which the radial wave functions are normalised to unity in the 
atomic sphere of radius S, i.e., f Q r 2 Rf(r)dr=l. In ASA, M i i+1 is expressed in terms 
of logarithmic derivatives Di=rR[/Ri evaluated at the sphere boundary. Skriver and 
Mertig derive the expression for M i)i+1 as 

M M+1 = - ( j) l (E F ) ( p l+1 (E F )x 

{[A - 0] [A+i +1 + 2)] + [E F - V(S)\ S 2 } [6) 

where V(S) is the one-electron potential and <fii(E F ) the sphere boundary amplitude of 
the I partial wave evaluated at E F . 

Numerical estimate to the magnetic energy are carried out using the fixed-spin- 
moment method In the fixed-spin-moment method the total energy is obtained 
for a given magnetization M, i.e., by fixing the numbers of electrons with up and down 
spins. In this case, the Fermi energies in the up and down spin bands are not equal 
to each other because the equilibrium condition would not be satisfied for arbitrary 
M. At the equilibrium M two Fermi energies will coincide with each other. The total 
magnetic energy becomes minimum or maximum at this value of M. Note that the 
two approaches, i.e., the self-consistent, floating-spin-moment method as well as the 
fixed-spin moment-method are equivalent in the sense that for a given lattice constant 
the magnetic moment calculated by the standard floating-spin moment approach is 
the same as the magnetic moment for which the fixed-spin moment total energy has its 
minimum |37| . In practice, the floating-spin moment approach sometimes runs into some 
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convergence problem. From experience, to avoid such predicaments in convergence, one 
may carefully monitor the mixing of the initial and final charges during the iterations 
and increase the number of k— points. Thus, for a better resolution to determine the 
change in the total energy with respect to the input magnetization, the k- mesh had 
1771 k— points in the irreducible wedge of the cubic Brillouin zone. 

By the fixed-spin-moment method the difference AE(M) (=E(M)—E(0)) for given 
values of M is calculated. The calculated AE(M) is fitted to the phenomenological 
Landau equation of phase transition which is given as 

AE{M) = ]T ^-a 2n M 2n (4) 

n>0 Zn 

for n = 3, where the sign of the coefficient a 2n for n = 1 determines the nature of the 
magnetic ground state, i.e., a 2 >0 refers to a paramagnetic ground state while a 2 < 
refers to a ferromagnetic phase. We have applied the approach described above to the 
study of carbon vacancy in MgCNi 3 [13j and 3d transition-metal-MgCiVz3 alloys [14]. 

3. Results and Discussion 

3.1. Equation of state 

Both X-ray and neutron diffraction techniques unambiguously report MgCNi 3 and 
ZnCNi 3 as cubic perovskites with their lattice constants determined as 7.201 and 6.918 
a.u., respectively. Assuming an underlying rigid cubic lattice, with Mg(Zn) at cube 
corners, Ni at the faces and C at the octahedral interstitial site, the total energy 
minimization were carried out to determine the equation of state parameters. The 
total energies calculated, self-consistently, for six lattice constants close to equilibrium 
were fed as input to a third-order Birch-Murnaghan equation of state [HEIIHS]- Note that 
the Birch-Murnaghan equation is derived from the theory of finite strain, by considering 
an elastic isotropic medium under isothermal compression, with the assumption that the 
pressure-volume relation remains linear. Hence, in the optimization procedure we have 
restricted the choice close to the equilibrium. 
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ZnCNi 3 


MgCNi 3 


a eq (a.u.) 


7.2255 


7.3041 


7.0558 


7.1387 


B eq (Mbar) 


0.3886 


0.3479 


0.4656 


0.4188 


B eq 


4.4106 


4.5255 


4.3444 


4.7813 



Table 2. Comparison of the equation of state parameters of cubic perovskite ZnCNi^ 
with that of MgCNi using the KKR-ASA method as described in the text. 

Since the choice of the exchange correlation potential in the Kohn-Sham 
Hamiltonian has proven to be sensitive in the structural characterization, we have carried 
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out the total energy minimization for two different approximations, namely the LDA and 
GGA j32l The results are shown in TableEJ The GGA considerably overestimates 
the lattice constant for either alloys, when compared to the experimental values. For 
MgCNi 3 , using the LDA description of the exchange-correlation, the lattice constant 
was calculated as 7.139 a.u., with the bulk modulus and its pressure derivative as 0.42 
Mbar and 4.78 respectively. These values are consistent with the earlier first-principles 
reports [UEI- The underestimation in the lattice constant for MgCNi 3 , however when 
compared to the experiments, owe to the over-binding effects in the LDA, and is a well 
known problem. 

For ZnCNi 3 the equilibrium lattice constant calculated using LDA yielded the 
value as 7.056 a.u. which when compared to the recent X-ray diffraction results [2] 
of 6.918 a.u., was found to be an overestimation. However, the results of the present 
calculations are consistent with the works of Johannes and Pickett [B] who employed the 
FP-LAPW method. Note that the consistency of the ASA calculations with that of the 
full-potential counterparts owe to the inclusion of the muffin-tin correction [2HJ - The 
KKR-ASA calculations further finds the bulk modulus and its pressure derivative of 
ZnCNi 3 as 0.46 Mbar and 4.34, respectively. As mentioned above, the overestimation 
of lattice constant in LDA is not so common, which suggests that the samples subjected 
to the experiments may be sub-stoichiometric. This was also emphasized by Johannes 
and Pickett (3| following the crystal structure characterization of MgC y Ni 3 alloys [E|. 
In the latter, both experiments [HI [T2j and theoretical calculations [TH] have shown that 
the lattice constant decreases as the C content in the material decreases. 

To look for the changes in the equation of state parameters as a function of C 
content in MgC y Ni 3 and ZnC y Ni 3 alloys, total energy minimization was carried out. 
The variation is shown in FigHJ For both MgC y Ni 3 and ZnC y Ni 3 alloys, the lattice 
constant as well as the bulk modulus decrease as C at% decreases. For MgC y Ni 3 
observed trend is consistent with the earlier X-ray diffraction measurements. The rate 
of decrease in the lattice constant is estimated as 0.142 a.u./at%C while for ZnC y Ni 3 , 
the lattice constant was found to decrease at the rate of 0.189 a.u. per at% of C . Though 
the lattice constant and bulk modulus showed similar trend for either alloys, the change 
in the pressure derivative of the bulk modulus as a function of y characteristically 
differed. 

The pressure derivative of the bulk modulus measures the rate at which the material 
becomes incompressible with increasing pressure, and is sensitive to the softness of the 
equation of state. In the Debye approximation for isotropic solids, which assumes a 
uniform dependence of the lattice frequencies with volume, one may express the average 
phonon frequency u> as B' eq cx fj^?, where V is the equilibrium volume of the unit cell. 
Note that volume for the vacancy-rich alloys decreases with decreasing y, while B eq 
maps a different trend for MgC y Ni 3 and ZnC y Ni 3 alloys. Such a behaviour indicates 
that the properties associated with the MgC y Ni 3 lattice could be characteristically 
different from that of the ZnC y Ni 3 counterparts. Also, one may note that the phonon 
spectrum for MgCNi 3 reveals that certain C modes play a vital role in the materials 
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Figure 1. The variation in the equation of state parameters, equilibrium lattice 
constant a eq (a.u.) , the bulk modulus B eq (Mbar) and the pressure derivative of 
the bulk modulus, as a function of y in AlgC y Ni 3 (open circles) and ZnC y Ni 3 (filled 
squares) calculated using the KKR-ASA-CPA method as described in the text. 



superconducting properties in addition to those of the Ni modes jUJl l4*T|. 

Partial replacement of Zn for Mg in MgCNi 3 has shown that the transition 
temperature decreases [42J. The findings also conclude that the nature of pairing 
mechanism in MgCNi 3 is conventional [42j. To study the changes brought about 
by Zn substitutions in the Mg sub-lattice of MgCNi 3 , we have carried out KKR- 
ASA-CPA calculations for Mgi_ x Zn x Ni 3 alloys. In FigJ2l we show the variation of 
the equation of state parameters of Mgi_ x Zn x Ni 3 alloys. The decrease in the lattice 
constant is consistent with the previous experimental report. The bulk modulus as 
well as its pressure derivative increases as x increases in Mgi_ x Zn x Ni 3 alloys. This 
clearly indicates that the average phonon frequency gets modulated as Zn replaces Mg 
in MgCNi 3 . 

3.2. Electronic structure 

In FiglBl we show the total and sub-lattice resolved partial densities of states of 
MgCNi 3 and ZnCNi 3 calculated at their respective equilibrium lattice constants. The 
characteristic features of both MgCNi 3 and ZnCNi 3 appear more or less similar with 
an exception of a sharp peak in the energy range —0.6 < E < —0.4, characteristic of Zn 
3d states. Being deep below on the energy scale compared to the Fermi energy, which 
is zero on the scale shown, one may expect Zn d states to be localized and thus behave 
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Figure 2. The variation in the equation of state parameters, equilibrium lattice 
constant a eq (a.u.) , the bulk modulus B eq (Mbar) and the pressure derivative of the 
bulk modulus, as a function of x in Mgi- x Zn x CNi3 calculated using the KKR-ASA- 
CPA method as described in the text. 



atomically. While for MgCNi 3 a small peak characteristic of the Mg-Ni bonding also 
appears in this energy range, but is less pronounced. Furthermore, the states near 
Ep are predominantly Ni 3d in character in both alloys, with little admixture of the 
C 2p states. One may also find that the position of the Ni 3d derived singularity is 
slightly lower in the energy scale for ZnCNi 3 than for MgCNi 3 , which is consistent 
with the previous results. The N(Ep) and the contributions to it from the sub-lattices 
are compared in Table QJ 
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14.557 
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1.199 


12.341 


4.509 


0.097 
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Table 4. Comparison of the total N(Ef) and sub-lattice resolved density of states of 
ZnCNi 3 and MgCNi 3 expressed in units of states/Ry atom. 



The reported values of N(Ep) for MgCNi 3 are at variance with the existing 
reports HJI3IE1IZII121I111I13I1EII1ZI- It appears that the value is sensitive to the 
basic approximations made in each type of the electronic structure method, and also 
to the parameters like that of the choice of Wigner-Seitz radii, choice of the exchange- 
correlation potential and others. However, under similar approximations, it is clear 
that for ZnCNi 3 the N(Ep) reduces by 12% in comparison with MgCNi 3 . This is 
consistent with the earlier first-principles FP-LAPW calculations [H]. The reduction in 
N(Er?) mav be lareelv due to the smaller lattice constant of ZnCNii. in comparison 



A first-principles comparison of the electronic properties of MgC y Ni% and ZnC y Ni 3 alloys 9 



50 




Energy (Ry) 

Figure 3. Comparison of the total and sub-lattice resolved partial density of states 
of ZnCNis and MgCNi^, calculated at their respective equilibrium lattice constants. 
The vertical line through energy zero represents the alloy Fermi energy. In the inset 
we show a blow up of the total density of states near the Fermi energy. 



more significant in MgCNi% when compared to ZnCNi^. This is consistent with the 
fact that a larger fraction of the charge would be transferred back to the Mg sub-lattice, 
in MgCNi 3 in comparison with that of the Ni sub-lattice. 

The change in the N(Ep) as a function of lattice constant in Mgi_ x Zn x CNi 3 alloys 
is shown in FigGB One may find that N(E F ) decreases for all values of x, with respect 
to lattice constant. However, N(Ep) as a function of x, at the equilibrium lattice 
constant, was found to deviate a little, as is evident from FigJHl This clearly suggests 
that the electronic structure properties are mainly governed by the CNi G octahedra. 
The atoms occupying the cube corners i.e., Mg and Zn, however, play a non-trivial role 
in determining the structural properties. 

3. 3. Hopfield parameter 

The Hopfield parameter r\ has been regarded as a local "chemical" property of an atom 
in a crystal. It has been emphasized earlier that the most significant single parameter in 
understanding the Tc of a conventional superconductor is the Hopfield parameter [35J. 
For strong-coupling systems, the variation in r\ is more important than the variation of 
(a; 2 ) in changing T c . Softening (u 2 ) often does enhance T c , but a significant change 
in the magnitude of Tc depends largely on a significant change in the r\ value rather 
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Figure 4. Comparison of the change in the total density of states near Fermi energy 
of ZnCNi-i and MgCNiz for a range of lattice constant as indicated in the figure. 
The vertical line through energy zero represents the alloy Fermi energy. 



than a small change in the corresponding (lo 2 ). As a matter of fact, we look for the 
changes in the rj from the three sub-lattices of these perovskites as a function of lattice 
constant as well as y in MgC y Ni 3 and ZnC y Ni 3 alloys. Note that for MgCNi 3 , it 
has been reported that the superconducting transition temperature Tc increases upon 
application of external pressure |48( I49j. Besides, experiments remain controversial on 
the strength of the electron-phonon interaction in MgCNi 3 EH EH IH21- It has 
been suggested that MgCNi 3 may be a strongly-coupled superconductor, however, the 
magnitude of Tc being marginally reduced due to the paramagnon interactions [i| 15 1|. 

In FigHHwe show the changes in the rj of MgCNi 3 and ZnCNi 3 as a function of 
lattice constant. It is clear from FigHHthat the r\ c and rjm linearly increase as a function 
of decreasing volume in either alloys. If the change in the average phonon frequency 
remains small, then either of these alloys could enhance the transition temperature with 
respect to external pressure. For MgCNi 3 this view is consistent with the previous 
experimental results. Similar characteristic feature holds for the vacancy-rich disordered 
alloys, the variation of which is shown in FiglTTl and fl~2l 

To have an understanding in the variation of r/c, r) vac and rjm where r] vac can be 
considered as the local chemical property of the electrons in the empty sphere, we show 
in FigUHl the change in these parameters as a function of y in both MgC y Ni 3 and 
ZnC y Ni 3 alloys. One may find that the variation of rj remains similar for both the 
alloys function of decreasing C content. 
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Figure 5. Comparison of the change in the total density of states of ZnC y Ni^ (solid 
line) and MgC y Ni^ (dashed lines) alloys calculated at their equilibrium lattice constant 
with y as indicated, he vertical line through energy zero in each panel show represents 
the Fermi energy. 



3.4- Magnetic properties 

Total energies from both the self-consistent, spin polarized and spin unpolarized 
calculations remain degenerate for MgCNi 3 and ZnCNi 3 alloys at their equilibrium 
lattice constants. This unambiguously shows that the materials are non-magnetic in 
nature. However, having suggested that MgCNi 3 is on the verge of a ferromagnetic 
instability (U El EU EH El], and also that incipient magnetism in the form of spin- 
fluctuations reside in the material, we attempt to compare the magnetic properties of 
MgCNi 3 and ZnCNi 3 alloys using the fixed-spin moment approach of alloy theory [25J. 

Numerical calculations of magnetic energy AE(M) for MgCNi 3 and ZnCNi 3 are 
carried out at over a range of lattice constants. The calculated results of AE(M) in the 
fixed-spin-moment method are shown in FigJHJ The calculated AE(M) curves are fit to 
the form of a power series of M 2n up to n = 3, for the polynomial as mentioned above. 
The variations of the coefficients, 02 in units of — , in -3-, and a$ in -3- as a function 

fl*B fig l^g 

of lattice constant are shown in Fig[03 The propensity of magnetism can be inferred 
from the sign of the coefficient which is quadratic in M, i.e., 02- The coefficient ai is 
the measure of the curvature and is positive definite when the total energy minimum is 
at M = 0, thus referring to a paramagnetic ground state. In general, when 02 becomes 
negative, it infers that there would exist a minimum in the AE — M curve at a value 
other than M = referring to a ferromagnetic phase at that value of M. The higher- 
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Figure 6. Comparison of the change in the sub-lattice resolved C 2p partial density of 
states of ZnCyNiz (solid line) and MgC y Ni^ (dashed lines) alloys calculated at their 
equilibrium lattice constant with y as indicated. The vertical line through energy zero 
in each panel represents the Fermi energy. 

order coefficients 04 and Oq however are significant and they control the variation of 
AE with respect to M. For example, for larger values of M, 0,4 and successively a e 
would dominate, and if 04(0-6) tends to be negative it would show a dip in the AE — M 
variation pointing towards a magnetic transition at a higher value of M. This, in the 
first-principles characterization of the magnetic properties of a material would refer 
to a possibility of a metastable phase at relatively large values of external magnetic 
fields. However, it has to be noted that calculations for large values of M can result in 
ambiguous results. Hence, it is suggested to carry out calculations for smaller values of 
M and use the above mentioned polynomial function up to the minimum order, where 
the curve fits with sufficient accuracy. 

Fig. Uni shows that for smaller values of lattice constant, the alloys show an 
enhanced paramagnetic character. One may also note that the variation in 04 and 
a§ coefficients are oppositely complimented and hence in the renormalized approach to 
include corrections due to spin-fluctuations, as suggested by Yamada and Terao [55J, 
they would cancel out in proportion preserving the trend in the variation of a?,. Thus, 
it becomes likely that the incipient magnetic properties associated with MgCNi 3 and 
ZnCNi 3 would decrease as a function of decreasing lattice constant. 
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Figure 7. Comparison of the change in the sub-lattice resolved Ni 3d partial density 
of states of ZnCyNij, and MgC y Nis alloys, over a small energy window around Fermi 
energy, calculated at their equilibrium lattice constant with y as indicated. The vertical 
line through energy zero in each panel represents the Fermi energy. 



4. Conclusions 

First-principles syudy of the electronic properties of MgCNi 3 and ZnCNi 3 , and also 

their non-stoichiometric alloys are carried out using the density-functional-based KKR- 

ASA method. We find that the lattice constant for ZnCNi 3 is overestimated, while for 

MgCNi 3 it is underestimated. This suggests that the material ZnCNi 3 subjected 

to experiments may be non-stoichiometric. As a function of decreasing C content 

in MgC y Ni 3 and ZnC y Ni 3 alloys, one finds an opposite trend in the variation of 

pressure derivative of the bulk modulus, which is proportional to the averaged phonon 

frequency. With electronic structure remaining essentially the same for MgC y Ni 3 and 

ZnC y Ni 3 , the results hint that non-stoichiometry may have opposite effects. Note that 

for 0.9 < y < 1.0, MgC y Ni 3 alloys are feebly superconducting, while according to 

the conjecture that has been made ZnC y Ni 3 is not. It can thus be inferred that the 

associated phonon modes in ZnCNi 3 and its disordered alloys may be characteristically 

different when compared to the MgCNi 3 counterparts. A comparison of the phonon 

spectra of these alloys thus become quite necessary to understand the absence of 

superconductivity in ZnCNi 3 , although it is iso-structural and iso-valent with MgCNi 3 . 
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